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ABSTRACT 


We present the first simulations of tidal stirring of dwarf galaxies in the Local 
Group carried out in a fully cosmological context. We use the ErisDARK cos¬ 
mological simulation of a MW-sized galaxy to identify some of the most massive 
subhalos > 10®Mq) that fall into the main host before z = 2. Subhalos are 

replaced before infall with high-resolution models of dwarf galaxies comprising a 
faint stellar disk embedded in a dark matter halo. The set of models contains 
cuspy halos as well as halos with ’’cored” prohles (with asymptotic inner slope 
7 = 0.6) consistent with recent results of hydrodynamical simulations of dwarf 
galaxy formation. The simulations are then run to = 0 with as many as 54 
million particles and resolution as small as ~ 4 pc using the new parallel N-Body 
code ChaNGa. The stellar components of all satellites are signihcantly affected 
by tidal stirring, losing stellar mass and undergoing a morphological transfor¬ 
mation towards a pressure supported spheroidal system. However, while some 
remnants with cuspy halos maintain signihcant rotational flattening and disk-like 
features, all the shallow halo models achieve v/a < 0.5 and round shapes typical 
of dSph satellites of the MW and M31. Mass loss is also enhanced in the latter, 
and remnants can reach luminosities and velocity dispersions as low as those of 
Ultra Faint Dwarfs (UFDs) In two cases a merger transforms the objects quickly 
into a spheroidal system already before infall. We argue that cuspy progenitors 
must be the exception rather than the rule among satellites of the MW since all 
the MW and M31 satellites in the luminosity range of our remnants are dSphs, 
a result matched only in the simulation with ’’cored” models. This suggests an 
early and efficient transformation of cusps into cores for Local Group satellites. 

Subject headings: Cosmological simulations — Galaxies: dwarf spheroidal — Galaxy 



evolution — Methods; numerical 
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1. Introduction 

Dwarf galaxy satellites of the Milky Way (MW) and M31 are the main focns of 
’’near-held cosmology”, providing a plethora of nsefnl constraints and tests for galaxy 
formation theories as well as for the natnre of dark matter in the ACDM paradigm. Yet 
their origin is still not completely nnderstood. Nearly all the known satellites, located 
within the expected virial radins of the bright spirals, i.e at i? < 300 kpc, are classihed 
as dwarf spheroidal galaxies (dSphs). They have fairly ronnd stellar components with low 


or negligible rotation bnt relatively high velocity dispersions (Mateo (1998); McConnachie 


(2012)). At the same time, nearly all of the dwarf galaxies located ontside snch radins and 


far from M31 are gas-rich dwarf irregnlar galaxies (dirrs), which have ongoing or recent star 
formation, battened disky shapes, and exhibit snbstantial rotational snpport {vrot/cy > 1). 
The different spatial distribntion of dwarf galaxy types is known as the morphology-density 
relation, and rehects an analogons effect seen in galaxy clusters and rich groups on larger 


mass scales (Dressier (1980)). 


Environmental mechanisms have been often invoked to transform dirrs into dSphs 


(Mateo (1998); McConnachie (2012)). One of the mechanisms proposed to achieve such 


transformation is tidal stirring, namely the cumulative effect of repeated tidal shocks 


occurring near the pericenter of the eccentric orbits of satellites (Mayer (2010)). 


The process has been investigated thoroughly with collisionless N-Body simulations 
adopting high resolution models of dwarf galaxies interacting with rigid or symmetric 


models of hosts on orbits with high eccentricities as those of cosmological subhalos (Mayer 


et ah (2001a); Mayer et ah (2001b) :Klimentowski et ah (2007); Kazantzidis et ah (2011a)). 


It has also been investigated in combination with ram pressure mass removal using 
hydrodynamical simulations of dwarf-main host interaction in which a diffuse gaseous halo 


consistent with observational constraints is added to the dark halo of the main host (Mayer 





































et al. (2006); Mayer (2010)) 


These studies established tidal stirring as a major mode of morphological evolution 
of dwarf galaxies, turning them from faint rotationally supported disks into even fainter 
pressure supported spheroidals, thus providing a natural explanation for the morphology- 
density relation ubiquitously observed in galaxy groups, including our own Local Group 


(Grebel et al. (1999)). 


Direct tidal heating as well as indirect heating stemming from bar and buckling 
instabilities triggered by tidal shocks have been both shown to play a role in tidal stirring, 
although the latter requires initial disks with mass-to-light ratios at the upper end of those 


observed in dwarf galaxies (Mayer et al. (2001b); Kazantzidis et al. (2011a)). These works 
also determined that ram pressure mass removal aided by tidal shocks can rapidly remove 
the gas content of satellites explaining why dSphs are devoid of gas and why they can have 
very high mass-to-light ratios. 

Recently the effect of the dark matter density prohle of dwarf galaxies on tidal 
stirring has been studied, motivated by the new hndings of cosmological hydrodynamical 
simulations of held dwarf galaxy formation which have shown that feedback from supernovae 
and massive stars can turn the cuspy halos predicted by AGDM into halos with shallower 
inner density slopes (7 in [0.3, 0.7] range), depending on the star formation rate achieved 


in the galaxy (Governato et al. (2010); Pontzen & Governato (2012); Shen et al. (2014); 


Governato et al. (2012); Di Gintio et al. (2013); Trujillo-Gomez et al. (2015)). There have 


also been attempts to show the existence of this effect specihcally for dwarf galaxy satellites, 
albeit in this case resolution is a potential issue ((Brooks et al. 2013)). [Lokas et al. (2012) 


and Kazantzidis et al. (2013) have thus updated the previous tidal stirring simulations by 


carrying out N-Body simulations in which disky dwarf galaxies have shallower dark matter 
prohles as those arising in the cosmological hydro simulations, and have checked results 
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against varying resolution. They have found that the transformation by tidal stirring is 
enhanced signihcantly when a shallower density prohle is adopted (7 = 0.2 or 0.6). They 
have also found that, since mass loss is also enhanced, remnants of dwarfs with shallow 


dark matter profiles can evolve into very faint objects such as Ultra Faint Dwarfs (Lokas 


et ah (2012)). The enhanced mass loss and transformation follow from the lower binding 
energy and more impulsive dynamical response of the system in the case of objects with a 
shallower halo prohle. 

However, these recent simulations, exactly like those of the last decade, are fairly 
idealized since they completely lack the cosmological context. Hence orbits are also 
idealized, being disconnected with the epoch of fall inside the host and with the shape of 
its gravitational potential, which is just assumed to be spherically symmetric. Instead, 
cosmological simulations show that the host is triaxial and grows substantially in mass and 


radius over the last ten billion years (Diemand et al. (2007)), which implies satellites can be 


accreted on increasingly wider orbits over time. Furthermore, in the idealized simulations 
subhalos can only interact with the main host, and not between them. Also, for simplicity 
these simulations have always explored a fixed or limited mass range for the satellites, rather 
than considering the wide range of masses from those of actual subhalos. It is expected 
that, the initial mass, and also the stellar mass relative to that of the halo, ought to have an 
important effect on the response of the system to tidal shocks. Furthermore, other studies 
have shown that mergers or strong interactions, while rare, can occur among dwarf galaxies 


just before infall, turning them into dwarf spheroidals very quickly (Kazantzidis et al. 


( 2011 b)), but leaving it completely open to identify the distinct signature of a remnant 
of one of such mergers as opposed to one from tidal stirring. Ideally, all these difficulties 
should be overcome in a fully cosmological hydro simulation of a MW sized galaxy, but with 
current resolutions around 100 pc and ~ 10 "^Mq in the baryons numerical results cannot be 
put on hrm grounds for objects that should be only a few hundreds parsecs in size. 
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Therefore, in order to make progress beyond the idealized simulation we devise a 
numerical strategy that in the past has only been attempted for simulations of galaxies 
in clusters ( Dubinski (1999); Mastropietro et al. ( 2005[ )). In this approach, known as the 
’’replacement technique”, cosmological subhalos are replaced with high resolution models of 
dwarf galaxies within the original cosmological volume. Models are constructed based on 
observational and theoretical constraints, and matched as closely as possible to the total 
virial mass of the parent subhalo. After the replacement the simulation is rerun to the 
present epoch. In this way galaxy satellites are automatically placed on realistic orbits 
and are subject to the full tidal interaction history of the parent subhalos. Due to the 
computational facilities and improved speed and parallelism of the codes, the simulations 
can span the equivalent of more than 10 Gyr of evolution with this technique, while in 
previous cluster simulations only a small time span, from z = 0.5 onwards, could be covered 


(Mastropietro et al. (2005)). In this paper we present the hrst replacement simulation 
carried out for a MW sized halo, using the ErisDARK cosmological simulation carried 


out in the concordance ACDM model (Kuhlen et al. (2013); Pillepich et al. (2014)). 

The simulations are carried out with state-of-the-art massively parallel N-Body codes, 
PKDGRAV2 (Stadel ( 2001| )) and GhaNGa ( Jetley et al.|| (2008); Jetley et al. (2010); Menon 


et al. (2015); http://www-hpcc.astro.washington.edu/tools/changa.html) on the GRAYXE6 


’Monte Rosa” of the Swiss National Supercomputing Gentre. 


The paper continues with Section 2 where the replacement technique and numerical 
simulations are presented in details. Next the analysis methods and results are enumerated 
in Section 3. Section 4 contains an extensive discussion regarding the implications of the 
results for the origin of dSphs, including their relation to observational data. Finally Section 
5 includes the conclusion and outlook of the presented work. 












































2. Numerical simulations and replacement technique 


2.1. Initial conditions and selection procedure 


The initial conditions of the simulations of this paper are based on the ErisDark run 


(Pillepich et ah (2014)), which is a dark matter only cosmological gravity run of 40 million 
particles using the ACDM cosmological model {Hq = 73 km s~^Mpc~^, Ug = 0.96,erg = 0.76, 
VLm = 0.268 and Ha = 0.732). The cosmological volume contains a 1 Mpe scale zoom-in 
region with three levels of rehnement that lead to a minimum mass of the dark matter 


particle of 1.2 ■ IO^Mq and a minimum softening length of 124 pc (Pillepich et al. (2014)). 
In the zoom-in region an isolated MW-sized halo assembles, with mass consistent with 
the lower limit for the estimated halo mass of the Milky Way, which provides a good £t 


with current kinematical data on the Milky Way stellar halo (Rashkov et al. (2013)). The 


companion hydrodynamical simulations, Eris and ErisLE, have shown to produce realistic 
spiral galaxies hj z = 0, which reproduce nearly all the main properties of the disk and 


bulge of the Milky Way (Guedes et al. (2011); Mayer (2012); Guedes et al. (2013); Bird 


et al. (2013)) 


Satellite galaxies falling into the main host after the last major merger, which takes 
place at z ^ 3, are considered for replacement with hi-resolution dwarf galaxy models 
including stars. We analize the entire available set of outputs of the original simulation with 


the use of the Amiga Halo Finder (AHF) (Gill et al. (2004):Knollmann & Knebe (2009)) 
with a virial overdensity criterion that varies with redshift according to the cosmology 


(Bryan & Norman (1998)). In this paper, in particular, we analyze the earliest infalling 


population of satellites among which there are survivors at z = 0. In the next paper we will 
present larger populations of dwarf galaxies which are accreted at later epochs. 


We thus proceed in the following way. First of all, among subhaloes with bound masses 


































greater than at z = 0, the satellite with the earliest infall is identified. The adopted 

threshold mass ensnres that we are considering well resolved objects, with at least 10^ 
particles (presnmably their progenitors are even better resolved at higher redshift, before 
tidal mass loss begins), to minimize nnmerical effects on their internal dynamical evolntion 


(Moore et ah (1996)). The object enters the host’s virial radins just before z = 2, when it 
is still unrelaxed, being the result of a merger between two subhalos. We trace the satellite 
back to z = 2.87 when the two progenitor merging subhalos are still well separated, namely 
their virial volumes do not overlap. Therefore the epoch corresponding to z = 2.87 is set 
as the epoch of replacement of the original ErisDark’s infalling subhalos with their high 
resolution versions. 

Next we identify other subhalos that at this epoch are located between 1 and 4 Ry^r of the 
host. We restrict the sample to the most massive ones, having pre-infall masses around 
IO^Mq and above. Among these we reject 60% that cannot be properly identified and 
replaced because they are either too elongated (axis ratio less than 0.65, possibly tidally 
disturbed by companions) or have a substantial fraction of their mass in subhalos. Four 
objects remain after this rejection procedure, with masses in the range 0.99 — 4 x lO^M©. 
Note that at z ~ 2 this translates into satellites having masses in the range 80-300 times 
smaller than the host, which has a virial mass of 3.00 x IO^^Mq and a virial radius of 55.6 
kpc at the corresponding epoch (both relatively high since ErisDARK has a fast assembly 


history with no more major mergers after z = 3, Pillepich et ah (2014)). We remark that 
the lower end of the satellite-to-host mass ratio is still in a regime in which dynamical 


friction is effective in eroding the orbital energy over many Gyrs (Colpi et ah (1999)), while 
this is a negligible effect in published idealized simulations that assume satellites of similar 


mass but within a host having a fixed mass IO^^Mq ([Kazantzidis et ah] (|2013[)). Finally, 


we observe that none of the four objects has a mass above lO^Mgoi at z = 0 in the original 
ErisDark run. Therefore they are not part of the largest subhaloes at the hnal epoch, which 
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can well exceed this mass at low z (Paper II, Pillepich et al. (2014)). In total we thns 
replace 6 snbhalos (fignre , of which 2 are merging during infall and the other four are 
accreted individually. The halo properties are listed in table [Tj 


2.2. Model galaxy generation 


The replacement models are generated using the GalactICS code (Kuijken & Dub inski 


(1995), Widrow & Dubinski (2005), Widrow et al. (2008)), which generates multi-component 
equilibrium models of galaxies. Each model galaxy is comprised of two components: the 
dark matter halo and the stellar disk. It is assumed by construction that the initial model 
has a stellar disk since this stems from the assumption that disky dwarfs are the progenitors 
of dSphs, the core idea of tidal stirring. The code creates halos with negligible angular 
momentum and a spherically symmetric density profile as described by the relation: 


Phalaij") 


2 ^ '^cr^ 1 I (r -Tc 

S^TTy-^ 



( 1 ) 


In the previous equation the scale radius and velocity scale are denoted with and cXs, 
and 7 stands for the cusp coefficient. In order to restrain the model to hnite dimensions 
a complementary error function with cut-off radius Tc and width is used. The velocity 
distribution is computed by solving the Eddington equation for the prescribed density profile. 


For each removed halo two replacements are created: one with 7 = 1.0 and another 
with 7 = 0.6. The cut-off radius and width are set to be linear functions of the virial radius 
of the original object: 

















-11 - 


- ‘T’ IT 

r. ' ??7.r ' 11 


( 2 ) 


= 


y/^^vir 

~w~ 


( 3 ) 


The scale radius is determined indirectly by assuming that the concentration of an 


average halo of mass Myir at redshift are given by the empirical relations found by Klypin 


et al. (2011): 


c{Myir,z) = Co(z] 
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( 4 ) 


With the concentration and virial radius set, the scale length can be determined. 


r. = 


( 5 ) 


The two variants corresponding to shallow and steep dark matter density prohle have by 
construct the same scale radius. In order to have the same mass enclosed in the same virial 
volumes with different cusp coefficients the scale velocity is adjust accordingly. The values 
of the characteristic quantities associated to dark matter halos with the virial quantities 
presented in table are listed in table 

The previously mentioned code can generate a disk component of a galaxy according 
to the cylindrical spatial distribution: 
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Pstars {r,z) = poexp ( - — 


seek ( — 


( 6 ) 


and cylindrical distribution of the radial dispersion: 
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(7) 


The reader who is interested in more details on how the code constructs the galactic 


model is referred to the publication Widrow et ah (2008). 


There are free parameters that have to be set using educated guesses based on available 
constraints on galactic structure at low mass scales. Most importantly, the mass of the 
stellar disk for each object is 1.2 x 10“^ of its virial mass Such ratio is chosen to reflect 
the stellar-to-halo mass ratio, inferred from the abundance matching analysis 


extrapolated to lower masses (Behroozi et ah (2013)). Past simulations of disk galaxy 


formation presented in the Kaufmann et ah (2007) paper associate to objects with maximum 


circular velocities in the 24-53 km s“^range, scale height values in the 0.08-0.18 kpc range. 
Considering the previously mentioned results the values of the scale heights of all the objects 
is set to the minimum of 0.08 kpc. The choice of constant scale height is motivated by the 
expectation that at the masses of interest the disk thickness acquires minimum threshold 
value. The value is dictated by the balance between the thermodynamics of disk formation, 
such as the effect of the cosmic ionizing flux providing a minimum temperature floor of 
~ 10^ K, corresponding to a bulk velocity dispersion ~ 8 — 10 km/s, and the gravitational 
potential of the halo. For all the stellar disks introduced in the simulation the scale height 
was set to 0.08 kpc. 

With the scale height and stellar mass to total mass ratio assigned, the scale length 
remains to be determined. The square of the scale length is naturally constrained to be 
proportional to the virial mass of the respective object. An observational study provided 


by Swaters et ah (2009) hnds several dwarf galaxies with disk scale lengths as low as 0.33 


kpc. Considering the inferred galaxy masses for the latter sample and the values assumed 
in previous work on tidal stirring based on various constraints (Kazantzidis et ah (2011a|)), 
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the reference value of Rd = QAkpc was set for an object with a virial mass of 10^Mgoi- Using 
such reference value as a normalization the following relation is used for the calculation of 
the scale length as a function of galaxy mass: 


Rd = 0.4 


(- 




0.5 


kpc. 


For each pair of objects the parameters of the stellar component are listed in table 


( 8 ) 


From the original cosmological simulation, quantities such as the virial mass and 
virial radius are measured (table [^. The virial quantities are used for determining all the 
parameters required for the generation of the pairs of high resolution galaxies that will be 
used as replacements. Namely the following quantities are determined: (cTsAsAcAw) the 
dark matter halo quantities and [Mstars-, Rd) the stellar disk quantities. The remaining free 
parameters are the number of dark matter particles and star particles. For all the models 
created both numbers were set to one million. 


2.3. Stability of the galaxies 

Each of the twelve high resolution generated objects evolve in isolation for 1 Gyr prior 
to their introduction in the cosmological environment. The state of the objects at the end 
of the isolation run can be seen in hgure and The purpose of these simulations is 
to verify the stability of the models and allow the system to remove instabilities prior to 
replacement. Initial discreteness noise is unavoidable, and it is known to seed waves in cold 


disk-like systems which could lead to mass redistribution (eg Mayer (2011)). These often 


manifest themselves as transient spiral waves, although the method adopted by our initial 
condition generator combined with the low self-gravity of the light disks embedded in a 
massive halo should minimize amplihcation of the waves. Furthermore, two-body relaxation 
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between more massive dark matter particles and lighter stellar particles may also lead to 
some spnrions evolution (Mayer et al. 2001). All isolation runs were performed using the 


pkdgrav2 code (Stadel (2001)). The softening lengths of the particles were set according to 
the equation: 


^dark 0.06 

^star 0.02 


(- 


M„ 


1/3 


M, \ 1/3 


^stars 


4 ■ lO^M,,; 


kpc, 

kpc, 


( 9 ) 

( 10 ) 


And are listed in table [H 


The baryonic distributions of the constructed galaxy models with a cusp coefficient 
equal to unity are presented in hgures and The previously mentioned images correspond 
to the objects at the end of the isolation run and at the moment of replacement. Figure]^ 
shows the cylindrically averaged surface density at the beginning of the isolation run and 
at the end. Small fluctuations are visible as small amplitude waves, which cause departure 
from perfect axisymmetry. Despite the existence of the local density fluctuations for both 
versions there are no signihcant changes happening on scales comparable to the scale length 
during the isolation run. Moreover, the analysis of the global mass distribution including 
the dominant dark matter component as shown by the circular velocity curves in hgure 
along with the radial distribution of the cumulative mass present in hgure do not show 
any appreciable changes. We conclude that the stability of the models is satisfactory. 


2.4. Replacement method and numerical simulations 

At the end of their evolution in isolation the galactic models replace the original objects 
in four steps. All bound particles of the original objects are removed after measuring their 
centre of mass (COM) velocity and position from the cosmological box at redshift 2.87. 
The particles of the newly created pairs of galaxies have their residual COM velocities and 
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positions removed. Afterwards their individual positions and velocities are rotated with 


random matrices generated with the Arvo (1991). Finally the COM phase space coordinates 
are translated to match the COM phase space coordinates of the original objects that they 
replace. The resulting modihed versions of the cosmological volume containing ErisDark are 
used as initial conditions for the gravity calculation runs down to redshift 0. ChaNGa was 
used to carry out these new cosmological runs. The code ChaNGa uses a fast hierachical 
oct-tree gravity solver, and is based on the Gharm++ parallel programming infrastructure 


(Kale & Krishnan (1996)). It leverages the object-based virtualization and data-driven 
style of computation inherent in Gharm-|--|- to adaptively overlap communication and 


computation and achieve high levels of resource utilization on large systems (Jetley et al. 


(2008), Jetley et al. ( 2010[ )). The Gharm-|--|- load balancing infrastructure is used to 
distribute pieces of this tree across processors. The leaf nodes are buckets that contain 
several particles (usually 8 to 32) whose force calculations are collectively optimized. At 
each level of the tree, multipoles are calculated to speed distant force evaluations. Time 
adaptivity is achieved by assigning individual timesteps to particles. A special load balancer 
was developed by the UIUG Gharm-|--|- team to efficiently handle the resulting processor 


load imbalance (Menon et al. (2015)). 


3. Results and data analysis 

The resulting raw simulation data is analized using the AHF in the same manner as 
performed with the original ErisDark simulation data. From the set of halos generated 
by the AHF code, the six pairs of satellites are identihed through the indices of bound 
star particles. As output of the AHF code, we use in our satellite analysis only the 
list of bound dark matter and star particles. There are cases when the halo hnder 
does not detect one of the objects at a particular epoch either due to the pericenter 




















shock or temporary interference with other large snbhalos. Snch events will appear as 
nnavoidable interrnption of the flow of data points on the analysis plots. As soon as the ob¬ 
ject is identified again by the AHF code, the flow of data points in the analysis plots resnmes. 


The center of the satellite is defined as being the position of the stellar density peak 
of the respective object. From this point a 3D half light/mass radius can be defined as the 
radius of the sphere centered on the previously mentioned density peak which includes and 
excludes half of the mass amount of bound stellar particles. All the values of the half light 
radius associated with each object at the available evaluation epochs are presented in figure 
The next step in the preparation of the satellite data for analysis is the alignment of the 
object. This is performed by defining a set of particles that include both dark matter and 
stars inside a sphere centered on the stellar density peak and with a radius equal to two 
half light radii. All the bound particles of the object including the ones exterior to the two 
half light radii sphere are rotated such that the moment of inertia tensor, corresponding to 
the set of stellar particles interior to the previously mentioned sphere, is diagonalized. This 
process is repeated for each individual object at each epoch of measurement. 

Once the satellite is identified and analysed there are a number of diagnostics that are 
routinely measured, as described in the following subsections. In addition to dynamical 
masses, dark matter masses and stellar masses (see tables and for the results at the 
starting and hnal time), the key diagnostics used to quantify the morphology of the model 
galaxies are the shape as measured by the axis ratios, the ratio between ID velocity 
dispersion a* and rotational velocity Vrot of the stars (note that we refer here to the actual 
stellar rotational velocity rather than to the circular velocity), and the surface density 
profiles of the stars. The characteristic velocities are measured inside (cr*) and around {vrot) 
the half-mass radius, which is defined for the whole 3D mass distribution. 
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Our classification scheme to decide whether or not the hnal remnant is a dSph is more 


stringent than the one used in (Kazantzidis et ah (2011a):Kazantzidis et ah (2013)). There 


two criteria were used simultaneously, namely that the object had to have c/a > 0.5 and 
Vrot/o'* < 1. With the latter criterion outliers from the typical structural properties of 
dSphs, such as the Tucana dSph, which has Vrot/o'* ~ 1, would be included. Here instead we 
use only the criterion Vrot/o'* < 0.5, which is indeed the case for nearly all observed dSphs 
and that, as we will see, automatically guarantees also that c/a < 0.5 in our simulated 
remnants at z = 0. Finally, for brevity we will refer to ’’cored models” when we talk about 
the models with inner slope of the dark halo 7 = 0 . 6 , and to ’’cuspy models” when we 
discuss models with an NFW prohle. 


3.1. Dimensions along the principal axes 

With the object aligned and principal axes identihed three maps of the projected stellar 
mass corresponding to the three pairs of principal axes can be generated. Perpendicular to 
each of the vectors dehning the axes directions a square grid with size 4 i?i /2 • 4 i?i /2 and a 
million square cells is dehned. For each cell a projected mass is dehned that equals the mass 
of all bound star particles whose projection on the grid is included in the respective cell. 
Using the information in the projected mass maps, the dimensions of the object along 
the principal axes can be determined. In the process of dimensions calculations ellipses 
containing cells with different levels of mass density are being htted. The distribution of 
these ellipses around the center of the image allows for the calculation of the ratio of the 
dimensions. For each of the three planes {x,y), {y,z), {z,x) the ratio between the shortest 
and longest dimensions c and a is measured: - , - , - 

^ cix,y^ az,x 


After the previously presented analysis is performed the hrst available scalar indicator 
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of the transformation of the object can be determined: 


c 

- = min 
a 



( 11 ) 


The time evolntion of the previonsly dehned qnantity is presented in hgnre 


3.2. Velocity dispersion and tangential velocity 


With the directions of the principal axes determined, the velocity dispersion is 
compnted in the direction of each axis. For the previonsly mentioned qnantities the stars 
considered for the calculations reside in a spherical volume of a radius Vc of interest, 


centered on the satellite. The velocity dispersion of the system associated with radius Vc is: 



( 12 ) 


In order to compute the rotational velocity associated to the system, three cylindrical shell 
volumes are dehned. Each of the volumes is dehned with a central axis constrained to pass 
through the maximum density region of the object and to be parallel to a principal axis 
dimension. The inner and the outer radii of the volumes are 0.95 Vc and 1.05 r^. 

For all the stars bound to the object and included in the previously dehned volume, 
the component of their velocity tangential to any concentric circle centered on the axis of 
the cylinder is considered for the average. The resulting average is the velocity uT, where 
ihj) e {y,z), (z^). 

With the dispersion and rotation quantities determined, the ratio: 



max ( 



(13) 


can be used as another indicator of the degree of transformation from rotationally supported 
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satellite to a dispersion supported galaxy (|^. In figure]^ the decrease of the ratio due to 
the effect of the tides can be observed. 


3.3. Projected mass surface densities and cylindrically averaged surface mass 

densities 


In addition to the quantities determined based on the bound amount of particles for the 
respective dwarf galaxies, the projected surface mass densities (PSMD) and cylindrically 
averaged surface mass densities (CASMD) have been determined. Unlike with the previous 
quantities where the nearby unbound particles are excluded from the calculations, in the 
PSMD and CASMD case they are included. The plots present in £gures|2} pT| [T^and 
13 are generated by projecting along the principal axes of the bound satellite all stellar 
particles inside a spherical volume of 10 kpc radius centered on the stellar density maximum 
of the galaxy. In the case of the CASMD calculations, the depth of the projection was 
reduced in order to include the stars with distances from the plane of projection smaller 
than 2.5 kpc. The resulting plots are available in figures |4, M 


3.4. Evolution of satellites 


3.4-1- Galaxy A 


Originally the most massive of the studied objects, satellite A enters the host halo on 
an orbit which brings it very close to the ErisDark’s center. Soon after the first pericenter 
passage, owing to the effect of dynamical friction the dwarf galaxy settles on an eccentric 
orbit with pericenter distances smaller than 5 kpc and apocenter distances less than 50 kpc 
Before the end of the simulation the cuspy version of the object completes 15 orbits, while 
the cored version completes 14 orbits, as it can be seen from the top left panel of figure 
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In both cases the satellite is found at z = 0 close to its maximum distance from the host. 
The three dimensional half mass radius decreases from the original 1300 pc values to 220 
pc and 110 pc for the NFW version and nonNFW version respectively. The time evolution 
of the half-mass radius can be inspected in figure 


The tides reach the stellar disc and unbinds 10% of the stellar mass before the second 
pericenter is reached. 90% of the stellar mass of the cored dwarf is lost prior to completion 


of the fourth orbit (hgure 16). In the case of the counterpart with steep density prohle, 
seven orbits have to be completed before loosing the same amount of stars. 


Based on the projected mass distribution of the satellite, it can be observed that the 
difference in the halo prohle has no clear impact on the evolution of the shape of the satellite 
(hgures 10 and 0 - The ratios between the shortest and longest dimension corresponding 
to both versions reaches the value 0.5 during the third pericenter passage and remains 
above this value throughout the simulation (hgure . 


From a kinematic perspective, the ratios between the rotational velocity and the 
dispersion associated with the respective half-light radius decrease to the 0.5 level diherently 
for the two versions (hgure |^. Before the hrst orbit is completed the cored satellite is 
described by a c/a of 0.5. At the same time the ratio of the NFW variant decreases to 
values in the 1.0 - 1.7 range. Later during the third orbit, the value 0.5 is reached as well. It 
should be noted that the hnal size of the object with 7 = 0.6 is comparable to the softening 
length of its dark matter particles. Therefore the very hnal stages in the evolution of the 
object should be regarded with caution as the response of the object to tides is ahected by 


the limited resolution ((Moore et al. 1996[ )). 

The evolutionary tracks of the descriptive parameters c/a, Urot/o'*, and M* are available 
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in figures |8j and respectively. 


3.4-2. Galaxy B 


Unlike the previous dwarf galaxy, the second object has a significantly wider orbit. 
The apocenter distance varies around the 100 kpc value with a pericenter distance less than 
10 kpc. During the 11 Gyr timespan of the simulation the satellite completes six orbits, 


and at the epoch corresponding to z = 0 it is aproaching its 7th apocenter (figure 16). The 
value of half mass radius shows large fluctuations in the case of the variant with shallow 
dark matter density profile, with large variations when the object is in the proximity of the 
periapsis. The half light radius of the alternative presents smaller fluctuations. Each of the 
versions of dwarf galaxies start with a half light radius of 1160 pc and end with a value of 
1010 pc in the case of the cored dwarf and 900 pc in the case of the cuspy dwarf (£gure[^. 
The stellar discs of the two versions start to lose significant amounts of stars during the 
second pericenter passage. By the end of the simulations the version with the higher central 
density has lost half of the initial stellar mass. In the case of the other variant this level of 
loss is reached before the 4th apocenter epoch. Moreover at z=0 90% of the original stellar 


mass is lost (figure 16) 


The velocity ratio corresponding to the NEW prohle reaches the level of 1.5 at redshift 0. 
At the same epoch the cored dwarf reaches a value of 0.5 (£gure|^. 

Spatially, the two versions of the object, maintain similar ratios of the minor and major 
dimensions until epoch 8 Gyr. Afterwards in the case of the variant with a cusp coefficient 
of 0.6, the ratio grows faster. The divergence in the velocity ratio appears at the same 
time; namely after the epoch corresponding to the 4th pericenter. Finally at z =0 , c/a 
fluctuates around 0.8 for the non-NFW variant and around 0.5 for the NEW variant (figure 
1^. When studying the projected density maps it is evident that the cupsy object maintains 
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a two-armed spiral pattern while the cored dwarf approaches spherical symmetry (hgures 


10 and 12). 


3.4-3. Galaxy C 


The third most massive object is originally on a trajectory with a pericenter distance 
of 40 kpc and an apocenter distance of 110 kpc. Soon after one orbit is completed it decays 
to a trajectory closer to the host. The following orbits are characterised by pericenter 
distances in the range 10-20 kpc and apocenter distances in the 90-100 kpc range. As the 


previous satellite before the last epoch it completes six orbits (figure 15). 


Kinematically, the characteristics of the two versions diverge after epoch 7.5 Gyr 
which roughly corresponds to the third apocenter passage (hgurej^. The indicator of the 
spheroidal transformation c/a clearly reflects the large difference in the transformation level 
(figure]^. While the cuspy object maintains a ratio around 1.5, the cored object reaches 
values below 0.2. The shape indicators of the two variants diverge earlier during the third 
pericenter passage and eventually converge to 0.7 values at z=0, roughly the time of the 
last apocenter passage. 


A visual inspection of the density map discloses a more complex case (figures 10 


and 12). The shape of the stellar distribution in the inner regions resemble case by case, 
the distribution of the satellite B pair. At the same time the outer, low density and 
spatially more extended regions are very different. Particularly they are more isotropically 
distributed in the case of object C than in the case of object B for both versions of the 
respective object. 
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The stellar masses of the two versions of satellite C remain approximately constant 
nntil the 7 Gyr epoch. From this moment on the stellar mass decreases from 1.7 10 ®Mq to 
1.0 IO^Mq in the case of the cnspy dwarf, and to 4 IO^Mq in its cored connterpart. The 
original valne of the half light radins for both versions was 810 pc, while the hnal valnes 
are 680 pc and 650 pc in the cnspy and cored dwarf, respectively (hgnrej^. Flnctnations 
similar to ones observed with object B are present, althongh with smaller amplitndes. 
Moreover they start to be visible at later epochs than in the case of object B. 


3.4-4- Galaxy D 


The object D is bonnd to ErisDark on an orbit with a pericenter distance of abont 20 
kpc and an apocenter distance in the 70-80 kpc range. Dnring the simnlation time frame 
the dwarf galaxy passes throngh 7 pericenters, and at z = 0 it is approaching the 8th (hgnre 
15). The tides of the host galaxy start to signihcantly affect the stellar component jnst 
before 7 Gyr ago. From this epoch onwards the cnspy object loses 3 x IO^Mq mass. Dnring 


the same period its cored connterpart loses 8 IO^Mq (hgnre 16). Similarly to the previons 
cases the stellar component in the cored model begins to lose snbstantial mass when the 
onset of large hnctnations in the valnes of the c/a ratio is also observed, rehecting the more 
impnlsive response to tidal shocks relative to the cnspy case (hgnre [^. The hnctnations are 
particnlarly enhanced for the cored version. Glose to the present epoch c/a reaches ~ 0.7 
for both variants of the model. 

Similarly, the velocity ratio starts to diverge at the 7 Gyr epoch and converge to similar 
valnes of 0.5 at z = 0 (hgnre [^. Shapewise the hnal state of the two versions appear to be 


similar (hgnres 11 and 13). Moreover the evolntion of the half light radins is similar for 


both cases, showing small hnctnations and a small decrease of the valne from the original 
640 pc valnes to 570 pc and 580 pc for the cnspy and cored object, respectively (hgnre [^. 
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3.4-5. Galaxy EF 


The Dwarf galaxy EF, which is the result of a merger between two objects comparable 
in size, evolves on a wide orbit with the pericenter distance varying from 30 to 50 kpc and 


apocenter distances with values as high as 160 kpc (hgure 15). Since infall the satellite 
completes three orbits. After the merger, the stellar mass of the object remains virtually 


constant for both variants (hgure 16). The ratio Uroi/f* stabilizes around 0.8 for the object 
resulting from the merging of two cuspy dwarfs. For the other variant, Urot/o'* continues to 
decrease until it reaches values close to 0.1 at 13 Gyr epoch (hgure |^. In the case of the 
shape ratio c/a the object resulting from the merger of the cuspy subhalos has on average a 
higher value; after 12 Gyr the average c/a is around 0.8 for the cuspy pair and 0.7 for the 
cored pair (hgure . 


The original values of the half mass radius prior to the merging of the two objects 
were 420 pc for object E and 390 pc for object F. Immediately after the merger the radius 
of the resulting objects is 550 pc and during the remaining simulation time it increases to 
610 pc, with no signihcant diherences between the cuspy and core variants (hgure [^. This 
highlights a clear diherence with tidally stirred dwarfs, in which the half-light radius always 
decreases. 


4. Discussion 

In the current work the evolution of hve pairs of satellites entering a MW-sized halo 
has been presented. Each pair consisted of two objects diherentiated by a single parameter, 
namely the cusp coefficient of the dark matter halo. In each pair, each of the objects was 
composed of a stellar disk component and a dark matter halo. The virial mass of the 
objects spanned one order of magnitude from 3.3 10®Mq to 3.83 IO^Mq. The simulation 
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commences at z =2.87 when all dwarf galaxies are in the proximity of their future host 
halo, yet outside its virial volume. The original distances to the host center ranged from 75 
to 165 kpc. Since the last major merger experienced by the host halo takes place before z 
= 3 the described population of objects is representative of the earliest infall objects and 
thus the most tidally affected population. Galaxies with later infall epochs will be affected 
by the tides on average to a lesser extent. 


All of the objects with cored density prohles match our dehnition of dSphs before 
z=0. At the same time only two out of four objects with cuspy density prohles satisfy the 
criterion. The more efficient transformation of cored models is expected as their response 
to tides is more impulsive owing to their lower gravitational binding prohle for comparable 


mass (Kazantzidis et ah (2013)). In particular the binding energy decreases much faster 


towards the center allowing the tides to disturb the stars found closer to the center of the 


satellite (Lokas et ah ( 2012 )) 


Figures [T^ and show the time evolution of various structural properties of 

all models, highlighting how the cored models in general ht better the range of values that 


classical dSphs and UFDs span for such properties (McConnachie (2012)). We caution, 
though, that there is a lack of objects with small half light radii, which may reflect the 
relatively small sample of initial conditions adopted here (see next section). 


The half mass radius, velocity and shape ratios have signihcant fluctuations on orbital 
timescales. Their fluctuations appear to be in phase with orbital position and grow in 
amplitude as more stellar mass is unbound from the respective satellite. In general the 
objects with an original 7 = 0.6 show signihcantly greater amplitudes in the fluctuations 
than the ones with 7 = 1 . 0 , again reflecting the different binding energy prohles which 
place cored models in the more impulsive regime for the response to tides (which will be 
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maximal near pericenter). The minimum of the fluctuations in half mass radius is reached 
when the object is near its periapsis, while the maxima is reached a few hundred million 
years afterwards. 

In general the satellites with 7 = 0.6 are also transformed faster into dwarf spheroidals 
than their counterparts. The final dark matter to stellar mass ratio are significantly higher 
in the cuspy remnants (see table [^, yet the hnal values in cored remnants are still high 
enough to match those of classical dSphs and UFDs. Assuming a stellar mass-to-light ratio 
(A//L). ~ 3 — 4, appropriate for an old stellar population, (table would yield values 
in the range M/L = 90 — 1000 for the cored remnants. Instead velocity dispersions and 
dynamical masses of cuspy remnants indeed appear to exceed observed values in model 
B and D, because of the very high central dark matter density. Note this result is not 
coincidental rather it reflects the fact that, while the core remnants have radii of less than 1 
kpc, the initial dark matter-to-stellar mass ratio was very high even inside the central kpc 
in order to match the constraints from abundance matching with the chosen density profiles 
of dark matter. 


All the tidally stirred remnants are found at i? < 100 kpc from the host at z = 0. 
Within such a distance from the MW and M31 all the galaxies in the mass range of the 
remnants (stellar masses 10^ — lO^M©) are dSphs (or UFDs, which structurally appear 


similar to classical dSphs, see McConnachie (2012)). One is thus tempted to conclude 
that this is evidence of the fact that the most massive satellites falling early inside the 
host, which is the type of objects we are studying here, have predominantly shallow dark 
matter prohles. In other words, dwarf galaxies with the properties of the cuspy remnants 
that did not transform fully into dSphs are not observed in the inner halo of the MW and 
M31, which means cuspy progenitors must be rare (the only irregular galaxies in the inner 
halo, the LMC and SMC, are a few magnitudes more luminous than any of our remnants). 
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Note that the latter inference is well grounded even with the small sample in this paper 
since specihcally this sample contains the objects that have experienced the longest tidal 
evolution possible inside the main host (early infall population). We thus expect a larger 
number of such incompletely transformed objects with disk-like features in the much larger 
sample with varying infall time that will be presented in Paper II. 


Note that idealized simulations already hinted at such a more efficient transformation 
of the cored models, still they never found such a strong effect because the mass ratio 
adopted between the satellite and the main host was at least an order of magnitude higher 
than in this work. The reason is that here the host halo is relatively low in mass (closer 
to ~ IO^^Mq than to IO^^Mq when the most massive satellites fall in at z ~ 2, an effect 
that only a fully cosmological simulation can account for. The case for cored progenitors is 
suggested by the notion that dwarf galaxy satellites which had the highest star formation 
rates, as required to reach a stellar mass in the range of 10 ® — IO^Mq already at z > 2, 
are those in which feedback was most effective to turn cusps into cores ([Governato et al 


(2010); [Governato et ah (2012); Shen et al. (2014); Onorbe et al. (2015)). These are the 
objects in which the star formation (SF) efficiency, dehned as ratio between stellar mass 


assembled and dark halo mass, was particularly high (Madau et al. (2014)). A high star 
formation efficiency in classical dSphs versus similarly low mass held dirrs is indeed seen in 
SF histories based on hi-res color magnitude diagrams (eg |Skillman et ah (2014); Aparicio 
et al. in preparation). 


Interestingly, gas-poor dwarfs with disk features have been detected in clusters (Lisker 


et al. (2006)) and have been explained as harassed galaxies that have accreted onto the 


cluster core only recently, thus have not been able to transform fully into dwarf spheroidals. 
It is tempting to speculate that such a population would be present in clusters even if the 
progenitors have dark matter cores before infall since a much larger fraction of galaxies 





























accrete at 2 ; < 1 in clusters as opposed to 2 ; > 1 for satellites into galaxy-sized halos. This is 
strongly suggested by the evolution of some of our models, which also in the cored variant 
transform fully into dSphs only after nearly 10 billion years (eg see models B and D in 


Moreover in order for the velocity ratio of a tidally affected object to reach the 0.5 
threshold for a spheroidal it must lose more than 25% of its original stellar mass (hgure 0- 
The total fractional mass loss was 1/3 at a minimum (in the cuspy version of model satellite 
D). This occurs for both cuspy and core variants, and implies a prominent tidal stream has 
to form if the object is transformed into a dSph. Detecting tidal features or tails might be 
difficult though because, as our simulations show, tails and streams show up at a surface 
density that is between 10^ and 10^ times lower than the central surface density of the 


object (14, 21), which corresponds to 7-10 magnitudes lower relative to the central surface 
brightness. The lowest values apply to cuspy models in which stellar removal is less efficient 


(figure 14 and 21). This confirms earlier results of Mayer et al. (2002). The transition to a 
marginally bound tidal tail and fully unbound stream following the orbital path is marked 
by the flattening of the profile at large radii, typically occurring at 3-6 kpc, or 3-6 times 


the half-light radius (figures 14 and 21), confirming past work with a variety of simulation 


techniques and models (eg Johnston et al. (1999); Mayer et al. (2002)). This means a very 
wide field of view would be needed to unequivocally detect the tidal debris surrounding the 
dwarf in addition to very deep photometry. The actual shape of the stellar profile varies 


also with projection and over time (figure 14 ), but the basic facts just highlighted hold. 
We also notice that in some cases the change of slope of the surface density profile between 
the inner bound and outer unbound stellar component is more evident in the cuspy variant, 
probably because in the cored models tidal mass loss occurs all the way to the very center, 
leading to more continuous removal of stars as a function of radius. 
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The evolution of the observable structural parameters for the simulated objects shown 


in hgures 18, 19 and hgure 20 also highlights that cored models can evolve enough in 


stellar mass and radius to end up in the region populated by Ultra Faint Dwarfs (UFDs), 


conhrming a previous claim of Lokas et ah (2012) obtained with idealized simulations. In 


general the objects evolve from the proximity of the locus populated by dwarf irregulars 
towards the locus of ultra-faint dwarf spheroidals by passing through the classical dwarf 
spheroidal locus. Except for the satellite A pair all galaxies can be found in the locus of 
classical dwarf spheroidals and the transition region between the faintest irregular dwarf 
galaxies and spheroidals. Rather than being necessarily reionization fossils, as it has been 


claimed in the literature (Ricotti (2010)), some UFDs might be thus just more extreme 
outcomes of tidal stirring. This scenario has the attractive aspect of allowing to explain 
naturally why there is a substantial continuum between the properties of classical dSphs 
and those of UFDs. Such a population of UFDs would have cored prohles, while on the 
contrary a population of pristine UFDs formed by gas collapse in mini-halos would likely 
maintain its original cuspy prohle as feedback from star formation would have a negligible 


effect (Penarrubia et ah (2012)). 


Despite the low hnal stellar mass such a population would yield no tension with the 


proposed core formation mechanism via feedback ( Penarrubia et ah| (2012)) because it 
would be the remains of a population of much more massive dwarf galaxies that were 
severely affected by the tides of the host. Indeed among our remnants the UFD-like object 
originates from one of the most massive infalling dwarfs, whose initial luminosity would 
be comparable to Fornax, exactly the mass scale at which the effect of feedback should be 


important (Governato et ah (2012); Shen et ah (2014)). Finding UFDs with cored prohles 
is thus a potential test of our scenario. UFDs forming from an early infall population 
of satellites would also be consistent with the notion that these galaxies typically lack 


an intermediate or young population (Brown et al. (2012)) as star formation would be 
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truncated soon after infall as gas is rapidly removed by tides and ram pressure when the 
cosmic ionizing background is still at its peak amplitude (Mayer et ah (2007)). 

The smallest two objects replaced merge while orbiting the host galaxy at signihcantly 
larger pericenter and apocenter distances in comparison to the other satellites. The effect 
of the tides is minimal with insignihcant stellar mass loss and small effects in the shape and 
velocity ratios for the merged objects with 7 = 1.0. In the case of the 7 = 0.6 mergers 
the velocity ratio further decreases due to the gravitational interaction with the host. If 
one considers the condition Vrotl<^* < 0.5 as the dehnition of a dwarf spheroidal, then the 
merging pair with 7 = 0.6 clearly reaches this state at the end. At the same time the merger 
remnant of galaxies with 7 = 1.0 fails the velocity criterion, although the later appears 
shapewise to be closer to a spheroid than the former. The formation of a signihcantly round 
object albeit with rotation still dynamically important cannot be caused by tidal stirring 
only since in that case the transformation of the initial disky shape is correlated with the 
redistribution of angular momentum due to tidally induced instabilities and tidal heating. 


(Mayer et ah (2007 


Hence hnding a dwarf with photometric properties analogous to those of classical dSphs 
but with a signihcant rotational velocity might be the signature of dwarf-dwarf mergers. 
Indeed low luminosity ellipticals, whose kinematics is reproduced well by mergers of disks 
with intermediate masses, do retain signihcant rotation (Cappellari et al. ( 2011 [ )). One of 
such rare objects might indeed be the Tucana dSph, which is on a very wide trajectory 


around the primary and exhibits appreciable rotation (Fraternali et al. (2009)). 


Figure M also shows that the change of slope in the case of the merger remnant occurs 
much more gradually, with a clear battening not appearing before 6-7 kpc from the center. 
This is not surprising as in in this case the morphological evolution is not accompanied 
by signihcant tidal mass removal. The diherent shape of the prohle relative to tidally 
stirred remnants is very interesting because it hints at a further possible way to distinguish 
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between a dSph formed by tidal stirring and one formed via dwarf galaxies mergers. The 
presence of a swift change from one exponential to another shonld be the signatnre of a 
tidally affected stellar system. 


5. Conclusion and Caveats 


To the authors’ knowledge the current work addresses for the hrst time the importance 
of the dark matter density prohle for the dynamical evolution and morphological 
transformation of disky satellites into spheroidals in a cosmological context. Therefore 
the work bridges the gap between detailed simulations of dwarf-host interaction and fully 
cosmological hydrodynamical simulations of MW-sized galaxies which do not have at the 
moment the resolution to reliably address this subject. This is an important step forward 
since the assembly history of the host, the infall orbits of the satellites and the structure of 
the potential well of the host, which all concur to determine the effect of tidal shocks, are 
all self-consistently taken into account. As a by product, dwarf-dwarf interactions are also 
take into account, and we witnessed a particular example of an interaction which gave rise 
to a merger. 


An important overall conclusion of this work is that tidal stirring is conhrmed to be 
an effective mechanism for transforming disky held dwarfs into objects with properties 
resembling dwarf spheroidals and even Ultra Faint Dwarfs (UFDs). Most importantly, it 
now appears clear that the morphological transformation induced by tidal stirring is a 
natural and inevitable consequence of hierarchical structure formation. The transformation 
is enhanced when the satellites have shallow dark matter density prohles, as in our cored 
models, conhrming previous results based on idealized simulations of individual dwarf-host 


interaction (Kazantzidis et ah (2013)). In the latter case, both tidal mass loss of the 


stellar component as well as tidal heating and internal redistribution are augmented as the 
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response of the galaxy is impulsive down to its central regions. As a consequence remnants 
can become as faint and puny as UFDs, which does not happen among the set of cuspy 
models. 


We also found that major mergers of satellite galaxies with shallow density profiles 
can effectively evolve into spheroidals on wide orbits. Since the merger happens before 
the satellites are accreted by the host, when they have lower velocities,it naturally leads 
to a dwarf spheroidal orbiting far from the center of the host. Especially in the cuspy 
case the remnant retains signihcant rotation despite achieving a spheroidal shape, which 
would not happen in tidally stirred remnants This is a plausible scenario for the origin of 
distant dSphs such as Tucana and Cetus, and confirms a previous claim based on idealized 
simulations (|Kazantzidis et aT (2011a)). 


Another important conclusion of this work is the identification of two signatures that 
could separate dwarf spheroidals formed due to a merger on large orbits from the ones 
formed due to the strong tides generated by the host galaxy. The hrst signature is the 
presence of a tidal stellar stream with a total mass comparable to the associated satellite. 
The second signature can be found in the stellar density/surface brightness profile of the 
object. It should appear as a fast transition on lengths of a few hundreds parsecs from 
the steep central exponential to a signihcantly shallower one. The objects that are not 
the source of such a stream or do not have a fast departure from the central brightness 
exponential decay most probably did not become spheroidal due to the gravity of their 
current host. They were probably formed through the merger of comparable galaxies. 


There are of course important caveats in the simulations, which primarily stem from 
the fact that they are carried out in the dissipationless approximation, namely without 
treating hydrodynamics, star formation and feedback. One caveat is the inexistence of the 
baryonic disk of the primary. This would increase signihcantly the central density of the 
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host, within 10-20 kpc, enhancing tidal mass loss and tidal stirring for satellites with orbits 


having pericenters smaller than ~ 10 — 20 kpc D’Onghia et ah (2010). Given the orbits of 


the satellites analysed in this paper, this effect would play a role at least for some of them 
(eg model A and model B, which have the smallest pericenters). The neglected effect should 
enhance the tidal transformation, but may lead to complete tidal disruption of objects 
that would have otherwise ended up as UFDs, or produce UFDs out of cuspy progenitors 


(Penarrubia et ah (2010)). Another caveat is the lack of a a gaseous component in the 


disky dwarf progenitors, but previous work has shown that this would be readily removed 
by ram pressure in the galactic diffuse corona over 1-3 Gyr for the mass range considered 
here {Vmax < 45 km/s), especially at z > 1 , when the cosmic photoionizing background 
maintains the ISM warm and loosely bound, thus suppressing further star formation even 


for the gas that has not been removed yet (Mayer (2010)). 


Furthermore, our modelling approach relies on the assumption that the progenitors 
of dwarf spheroidal satellites started out their life as disky dwarf galaxies. Gosmological 
hydrodynamical simulations show that this is a well justified assumption for dwarf galaxies 


forming in halos in the range 10® — 10 ^®Mq ( Governato et al)| (2010); Shen et ah (2014)), 
while at lower masses field dwarfs can have a rather thick and turbulent disk, with a ratio 
l^rot/ rsj 1 (Shen et ah, in preparation), which appears consistent with the kinematics 


of the faintest dirrs in the Local Group (Mateo (1998); McGonnachie (2012)). A thicker. 


kinematically hotter disk could have a different response to tidal shocks. On the one hand 
the triggering of non-axisymmetric instabilities, such as bars and spiral arms, would be 
weaker due to the reduced self-gravity, on the other hand the disk binding energy will be 


lower, so that the direct effect of tidal heating could be stronger. Kazantzidis et al. (2011a) 


have considered thicker disks in the wide range of initial conditions they explored in their 
tidal stirring simulations, hnding small differences in the transformation as the initial disk 
thickness was varied by a factor 3, with only a slight tendency of thicker disks to transition 
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faster to a nearly isotropic configuration. However, future cosmological hydrodynamical 
simulations with enough resolution to study the morphological evolution of satellites with 
the same level of detail as done here will be required to assess the impact of disk thickness 
and kinematics, of both gas and stars, on the efficiency of the transformation. 

Finally, here we have considered a relatively small sample of satellites as we focused 
on the earliest infalling objects with a sizable mass, large enough to host luminous dwarfs 
and be sufficiently well resolved. The advantage of the small sample is that we could study 
in detail the different dynamical evolution between cuspy and cored variants of the same 
object. In addition, since these early infallers are exposed for the maximum amount of 
time possible to the tides of the primary we can safely assume that the results provide an 
upper limit on how deeply the structure of the dwarfs can be modihed by tidal shocks in 
cored versus cuspy variants. As a result, we inferred that the population of satellites of 
the Milky Way should be dominated by remnants of cored objects otherwise we would see 
a population of dwarfs with intermediate properties between dwarf irregulars and dwarf 
spheroidals at ~ 100 — 200 kpc galactocentric distances, which is not seen neither in the 
MW or the M31 halo. In a forthcoming paper we will study a much larger sample of 
satellites which includes objects that are accreted at low redshifts. The size of the sample 
that we will present will be of the same order of the number of luminous satellites known 
so far to he within the putative virial radius of the Milky Way, thus allowing to make 
statistically meaningful statements on the evolution of the entire population of satellites. 

The first author is grateful to Annalisa Pillepich, Davide Fiacconi, Michela Mapelli, 
Doug Potter, Stelios Kazantzidis and Rok Roskar for their help. Funding for the work 
was provided by the Institutes for Theoretical Physics and Computational Science at the 
University of Zurich. 
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Table 1. The virial quantities of the replaced halos and their distance from the center of 

the host. 


Satellite 


Rvir [^_pc] 

dhost[kpc] 

A 

3.83 

13.0 

140 

B 

3.04 

12.0 

75 

C 

1.51 

9.5 

82 

D 

0.94 

8.1 

157 

E 

0.40 

6.1 

165 

F 

0.33 

5.7 

159 
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Table 2. Scale quantities and concentrations of the different halos and the properties of 

their associated stellar disks. 


Satellite 

A 

B 

C 

D 

E 

F 

c 

6.0 

6.1 

6.3 

6.4 

6.8 

6.8 

rs [kpc] 

2.18 

2.00 

1.53 

1.27 

0.91 

0.85 


88.4 

82.0 

65.2 

55.9 

42.4 

40.0 


84.7 

78.6 

62.4 

53.4 

40.4 

38.0 

Mgtars [IO^Mq] 

4.50 

3.57 

1.75 

1.08 

0.44 

0.37 

Rd\pc\ 

80 

80 

80 

80 

80 

80 

hd\pc] 

10 

10 

10 

10 

10 

10 
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Table 3. Masses and force resolutions of dark matter and stellar particles. 


Satellite 

Msp [Msoi] 

Mdmp [Afso/] 

^sp\pc] 

(^dmp [pc] 

A 

4.50 

3830 

9.6 

93.8 

B 

3.57 

3040 

8.9 

86.9 

C 

1.75 

1510 

7.1 

68.8 

D 

1.08 

940 

6.0 

58.8 

E 

0.44 

400 

4.5 

44.2 

F 

0.37 

330 

4.2 

41.6 
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Table 4. Bound stellar and total masses inside the 500 pc radius spheres centered on the 
density peaks of the satellite galaxies with steep (left) and shallow (right) density prohles 

at z = 2.87. 


Satellite 

"^50o[10'^Mq] 

"^500 [10^ 

"^50o[10’^Mq] 


A 

8.75 

6.26 

4.32 

6.09 

B 

7.98 

5.91 

3.99 

5.82 

C 

6.08 

4.88 

3.37 

4.80 

D 

4.96 

4.06 

2.91 

4.09 

E 

3.44 

2.67 

2.21 

2.69 

F 

3.17 

2.40 

2.06 

2.40 



Table 5. Characteristic masses and mass ratios at last epoch for the galaxies with 
original 7 = 1.0 (npper rows) and 7 = 0.6 (lower rows) 


Satellite 

A 

B 

C 

D 

EF 

^sooIIO’^Mq] 

0.64 

3.62 

3.21 

3.07 

4.44 

^500 [10^ "^ 0 ] 

0.25 

3.49 

3.23 

3.37 

3.50 


256 

104 

99 

91 

127 

^h^^und/ ^hound 

238 

102 

112 

128 

370 

^soollO’^M©] 

0.14 

0.39 

0.60 

0.71 

2.61 

^500 [10^^©] 

0.06 

0.72 

1.42 

1.75 

3.42 


234 

54 

42 

41 

76 


234 

39 

42 

47 

312 
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Table 6. Characteristic velocities and velocity ratio at the last epoch for the galaxies 
with original 7 = 1.0 (upper rows) and 7 = 0.6 (lower rows) 


Satellite 

A 

B 

C 

D 

EF 

(T*[/cm/s] 

6.20 

11.04 

9.88 

10.01 

11.46 

Vcirc[km/s\ 

9.90 

18.85 

17.33 

16.99 

20.96 

^-k / ^circ 

0.63 

0.59 

0.57 

0.59 

0.55 

ai,[km/s\ 

- 

3.52 

4.54 

4.92 

8.95 

Vcirc[km/s\ 

- 

6.0 

8.0 

8.1 

17.8 

/ ^circ 

- 

0.59 

0.57 

0.61 

0.50 
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X [cMpc] 

Fig. 1.— The projected density map of the ErisDark cosmological simulation at redshift 
2.87 for a cubic comoving megaparsec section of the zoom-in region with the same depth is 
presented. The replaced halos are visible due to the presence of the stars marked with red 
to yellow colors in the image. 



[kpc] X [kpc] X [kpc] 


42 



y [kpc] z [kpc] X [kpc] 



log ) 


Fig. 2.— The collage of images describes the stellar distribution of the high resolution 
version of satellites A, B, and C prior to the replacement. 
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Fig. 3.— The collage of images describes the stellar distribution of the high resolution 
version of satellites D, E, and F prior to the replacement. 
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Fig. 4.— Cylindrically averaged stellar surface densities of the high resolution objects before 
replacement. In the hrst/second row the halos with 7 = 1.0/0 .6 are presented. In the left 
column the objects before the isolation run, while in the right the objects after the isolation 
runs are presented. The steeper drop in the value of the function for low stellar densities 
is caused by the cuttoff of the stellar disk which is initiated at hve scale lengths from the 


center. 
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Fig. 5.— Spherically averaged circular velocities of the high resolution objects as a function 
of radius. In the hrst/second row the halos with 7 = 1 . 0 / 0.6 are presented. In the left 
column the curves corresponding to the objects before the isolation run are shown, while in 
the right the curves corresponding to the objects after the isolation runs are presented. The 
solid color vertical lines mark the radii corresponding to the maximum circular velocities 
and the dashed lines mark the force resolutions of the largest and smallest objects. 
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Fig. 6.— In the left row the cumulative masses of the generated models prior to the isolation 
run are presented as a function of distance from the center before the isolation run. In the 
right column the values of same functions after the isolation run are presented. The top row 
corresponds to the objects with 7 = 1.0 while the bottom to the objects with 7 = 0.6. The 
colored vertical lines correspond to the original virial radii of the respective halos. 
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Fig. 7.— The collage presents the time evolution of the 3D half light/mass radii for the five 
resulting pairs of satellites. The data points corresponding to galaxies with shallow density 
prohles are marked in blue and the points corresponding to galaxies with steep density 
prohles are marked in red. 
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Fig. 8.— The evolution of the dimension ratio c/a of the objects during the simulation. 
Red/blue dots represent the satellites with steep/shallow central density profile. The gray 
region of the last subplot marks the epochs when the merging process is in progress according 
to the half light radius indicator of|7| 
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Fig. 9.— The evolution of the ratio between the circular velocity and the velocity dispersion, 
both estimated for the 3D half-light radius.The data points corresponding to galaxies with 
shallow density prohles are marked in blue and the points corresponding to galaxies with 
steep density prohles are marked in red. 
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Fig. 10.— The projected surface mass density maps of objects A (first row), B (second row), 
and C (last row) with 7 = 1.0 at redshift 0. The first/second/third column corresponds to 
the projections along the smallest/middle/largest dimension. 
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Fig. 11.— The projected surface mass density maps of the objects D (hrst row) and EF 
with 7 = 1.0 at redshift 0. The hrst/second/third column corresponds to the projections 
along the smallest/middle/largest dimension. 
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Fig. 12.— The projected surface mass density maps of objects A (first row), B (second row), 
and C (last row) with 7 = 0.6 at redshift 0. The first/second/third column corresponds to 
the projections along the smallest/middle/largest dimension. 
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Fig. 13.— The projected surface mass density maps of the objects D (hrst row) and EF 
with 7 = 0.6 at redshift 0. The hrst/second/third column corresponds to the projections 
along the smallest/middle/largest dimension. 
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Fig. 14.— The cylindrically averaged surface mass density profiles are presented. All 
quantities are computed for the direction corresponding to the shortest dimension. Red/blue 
dots mark the satellites with steep/shallow central density prohle. The top/middle/bottom 
row illustrates the surface density of satellites A and B / C and D / EF. 
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Fig. 15.— Distances between the satellite galaxies and their host halo. The versions of 
dwarf galaxies with 7 = 1.0 / 0.6 are marked in red/blne. For object EF (the merging case) 
the gray region marks the epochs before the evolntion of the half mass radius stabilizes (as 


can be seen in figure . 
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Fig. 16.— The change in the mass of the bound stellar component over time is presented. 
Red/blue dots represent the satellites with steep/shallow central density prohle. 
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Fig. 17.— The correlations between the amonnt of stellar mass lost and velocity ratios. 
Red/blne dots represent the satellites with steep/shallow central density prohle. From top 
to bottom the snbplots correspond to satellite (A and B), (C and D) and EF. 
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Fig. 18.— The evolution of the mass of the bound stars and velocity dispersion of simulated 
satellites in comparison to the values corresponding to the local dwarf galaxies as compiled 


in the paper McConnachie (2012). Red/blue dots represent the satellites with steep/shallow 


central density prohle. From top to bottom the subplots correspond to satellite (A and B), 
(C and D) and EF. The black rhombuses mark the local dwarf spheroidals found within 300 
kpc from MW’s center. 
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Fig. 19.— The evolution of the mass of the bound stars and the half light/mass radius of 
the simulated satellites in comparison to the values corresponding to the local dwarf galaxies 


as compiled in the paper McConnachie (2012). Red/blue dots represent the satellites with 


steep/shallow central density prohle. The black rhombuses mark the local dwarf spheroidals 
found whithin 300 kpc from MW’s center. The green squares mark the held dwarf irregulars. 
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Fig. 20.— The evolution of the mass of the bound stars and dynamical mass within the half 
light radius of simulated satellites in comparison to the values corresponding to the local 


dwarf galaxies as compiled in the paper McConnachie (2012). Red/blue dots represent the 
satellites with steep/shallow central density prohle. The black rhombuses mark the local 
dwarf spheroidals found whithin 300 kpc from MW’s center. The green squares mark the 
held dwarf irregulars. 
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Fig. 21.— The cylindrically averaged surface mass density profiles are presented. The top 
row illustrates the surface density of satellite B computed along all three principal axes. The 
bottom row displays the surface density of satellite C at different moments in it’s orbital 
movement. Left column corresponds to objects with an original cusp coefficient equal to 1.0. 
Right column corresponds to objects with a cusp coefficient equal to 0.6. 
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